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ABSTRACT 


This investigation analyzed families of ascent 
trajectories to determine total velocity.require- 
ments, division of velocity impulses and the ex- 
pected errors in position and velocity upon arrival 
in the terminal region caused by uncertainties in 
the position, magnitude and direction of the 
velocity vector at burnout of the booster stages. 


Conditions at burnout of the booster and apogee 
of the trajectory were solved using an IBM 650 
computer; the error analysis was accomplished by 
numerical integration on IBM 709 equipment repeated 
up to the apogee time of the zero-error trajectory. 


From the tabulated and graphical results pre- 
sented herein, some trade-offs between propulsion 


requirements and the accuracy of guidance and con- 
trol systems can be Seen. 


Thesis Supervisor: Paul E. Sandorff 


Title: Professor of Aeronautics and Astronautics 
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OBJECT 


It is the object of tnis thesis to 
present in grapnical form the effects 
of errors in position and velocity at 
booster burnout on the position and 
velocity of a rendezvous venicle at 


zero-error apogee time. 





CHAPTER I 


INTRODUCTION 


Requirements 
Long before the launching of tne first earth 


satellite the military, scientific, ecenomic and 
political ramifications of such a vehicle were con- 
sidered and discussed. Each ef the above considera- 
tions generates adequate reasens for a satellite 
rendezvous. 

The most important and prebably the most 
immediate requirement is in the field of the mili- 
tary. In orbit about the earth today are satellites 
capable of reconnaissance; little additional equip- 
ment is required to extend this capability te ferret 
and jamming missions or ICBM early warning. Inter- 
rogation and possible negatien ef such enemy vehicles 
are missions that require & satellite rendezvous. 

Satellite stations in space require the physical 
exchange ef persennel, equipment and fuel. Com- 
pletion of such an exchange requires a "seft" rendezvous 
and centinuation of the twe vehicles in the same orbit 


until the transfer is complete. 





Beth civil and military erganizations are 
developing communications satellites; these de- 
vices invelve much complex and expensive equip- 
ment that vill, in time, require repair er re- 
placement. Satellite rendezveus is a solution 
to tnis problem. 

Communications satellites are given Special 
consideration in this report by the inclusion of 
data for orbit altitudes as great as 20,000 nauti- 
cal miles above the surface of the earth. Earth 
satellites much above this altitude are of little 
interest since they are subject to large pertur- 
bations by the sun and moon and hold no advantage 
over the satellite with a period of 24 hours. The 
effect of such perturbations and the difficulties 
to be encountered in maintaining a high altitude 


orbit are discussed by Edelbaum in reference 1. 


Background 

The Subject of satellite rendezvous has re- 
ceived its share of attention in the past several 
years. Literally hundreds of papers and articles 
have been written on the various phases of rendez- 
vous and interception. Most of these works are 


very general and deal primarily with ideas rather 
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than figures and hardware. Other authors look into 
the problem in detail and find it to be enormous; 
it seems to take a well-funded project to generate 
any amount of quantitative date on a particular 
phase. Many Studies are initiated on the basis 
that the proposed payload and vehicle character- 
istics are known and the task is to find out how 
(or if) the desired iis wen be accomplished. 
Much of this work is security classified and is, 
therefore, unaccessible. It is the purpose of 
this thesis to conduct an academic yet practical 
study of satellite rendezvous so that students may 
gain an insight into a very complex problem. 

An important part of this subject that is 
specifically not included in this thesis is the 
topic of guidance and control during the terminal 
phase of the rendezvous maneuver. This topic 
appears to be a study in itself inasmuch as the 
terminal phase shifts to "closed loop" control 
with feedback from radar as the rendezvous vehicle 
closes on the target. Some facts about the capa- 
bilities of the acquisition system and, in fact, 
entire specification of the means to be used by 
the rendezvous vehicle in completing the terminal 


phase must be outlined to justify the definition 





2 
of a successful arrival. Sears and Felleman in re- 
ference 2 give an excellent example of this type 
study. Nason establishes the conceptual limitation 


of proportional navigation in reference 5. 





CHAPTER II 


STATEMENT OP THE PROBLEM 


General 

The earth satellite ls a unique device 
since, unlike vehicles that operate in the earth's 
atmosphere, its flight is a regular pattern in 
Space determined primarily by the velocity imparted 
to the vehicle at the beginning of its orbit. While 
velocities of earth satellites are high compared to 
any motion within the earth's atmosphere, the motion 
of the satellite is predicted with accuracy and it 
is this statement that forms an important base uti- 
lized in this thesis. 

Satellite rendezvous is the problem of bringing 
together two objects of relatively small size when 
the original distance between the two objects and 
their velocities are both very large by ordinary 
engineering standards. When the target orbit ele- 
ments are known, the task of the rendezvous vehicle 


(RV) is to match them. In practically any circum- 





stances, even the simplest case, there exist an 
infinite number of possibilities for accomplish- 
ing a rendezvous. 

The general case of matching the ephemeris 
of an arbitrary target after launch from an 
arbitrary point in space (surface of the earth) 
is unmanageable and probably unrealistic. Al- 
though an effort is made to keep tne problem as 
general as possible, certain conditions are set 


before proceeding. 


Ine Specific Problem 
The study begins with the assumption that the 


orbit of the target and the trajectory of the RV 
are in the same plane. If the orbit of the target 
has an inclination less than the latitude of the 
launching point then plane-change is required as 
the target never passes over the launch point. 
With angles of inclination greater than this lati- 
tude the launch point is in the plane of the target 
twice each day giving two opportunities for launch 
into the target plane. Except for polar orbits, 
one of these opportunities requires a westward 
launch, an undesirable choice, thus the desirable 
eastwará launch into the target plane is available 


only once a day. 





Now consider further the odds that the target 
is at a suitable position for rendezvous at the | 
time the launch point passes through the target 
orbit. Obviously a wait of many days may be re- 
quired before the correct conditions are met, 

A study of the wait time problem is made by 
Houbolt in reference 4. In this paper he de- 
monstrates for a specific case the limited number 
of rendezvous opportunities available in a given 
period and points out the increase in this number 
given by the "Adjacent Launching Technique". In 
this scheme the RV is launched into an orbit of 
the same altitude as the target with the target 
at a suitable position for rendezvous. When the 
RV crosses the target orbit it adapts to this 
new orbit with a plane-change maneuver and begins 
the terminal phase of the rendezvous. 

inasmuch as changes of orbit inclination are 
accomplished by an impulse normal to the original 
velocity vector, the change of plane requirement 
is in addition to the in-plane requirements and 
is not considered further. 

Since the target is, by definition, the point 
in space to which we wish to move the RV, it is 


reasonable. to assume that the target has a regular 





motion easily predicted. A circular orbit best 
meets this requirement and this is not an un- 
realistic assumption; for many reasons an approxi- 
mately circular orbit is desired for reconnaissance, 
survey or communications satellites. 

Now within a given plane and with a circular 
target orbit there still remain an unlimited 
number of trajectories that meet the requirement 
of placing the RV at the satellite's position in 
Space. These trajectories form the basic conic 
sections, ellipse, parabola and hyperbola, passing 
through the two required points. Only a minimum 
of computation is required to see that hyperbolic 
and parabolic trajectories have excessive velocity 
requirements when compared with the elliptical 
path. The elliptical trajectories that qualify are 
those with apogee greater than or equal to the al- 
titude of the target. The former, often called a 
"fast ellipse", crosses the target orbit and this 


is undesirable for two reasons: 


1. To match the orbit of the target the RV 
requires in-plane normal thrusting at the cross- 
over point. Normal thrusting is inefficient and 


is only used as a last resort. 
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2. Closing velocities at the cross-over are 
very high and rendezvous requires critical timing. 
The ascent ellipses considered in this thesis 
are those with apogee equal to the altitude of the 
target. This type trajectory gives a flight path 
at apogee that is tangent to the orbit of the 
target. 
| Within the conditions of being in-plane and 
tangent at arrival, there is an alternative, four- 
impulse, "parking orbit" ascent. Briefly, parking 
orbit ascents consist of a launch into a low aiti- 
tude orbit, a waiting period until the relative 
positions of RV and target are correct, then a 
transfer to the target orbit. The primary advantage 
of such a scheme is tne relaxation of strict re- 
quirements for leaving the launching pad at a given 
time; the low altitude orbit allows the RV to "gain" 
on the target. The disadvantages of such a system 
are related to this thesis in that errors in burnout 
velocity are not only multiplied by the requirement 
for two "transfer" maneuvers, but initial errors in 
many cases "grow" during the wait period. These 
and other factors are considered in detail by Duke, 


Goldberg and Pfeffer in reference 5. 





11 

Coasting in free flight from burnout of the 
booster stages until the RV is at the altitude 
of the target is termed a "direct ascent". Such 
a trajectory has a minimum velocity requirement 
in the classical Hohmann maneuver and an upper 
limit in the vertical ascent. This thesis will 
be concerned with these two-impulse type rendez- 
vous trajectories. 

To review, then, the specific problem in- 
vestigated in this thesis is an in-plane rendezvous 
with a target in circular orbit utilizing a direct 
ascent elliptical trajectory that is tangent to the 
en orbit. 

Basic Assumptions 

The earth is assumed to be a homogeneous 
sphere. The earth is, of course, quite irregular; 
it is more an ellipsoid than a sphere. 

This oblateness has three main perturbating 
effects upon the orbit of an earth satellite, all 
dependent upon altitude and inclination of the 
orbit: 

1. Regression of the nodes. This is motion 
of the orbítal plane about the earth's axis in a 
direction opposite to that of the satellite. 





12 

2. Precession of the line of apsides. This 
motion is an angular rotation of the axis of the 
orbit in inertial space. It may be in the same 
or the opposite direction of the satellite motion 
depending upon tne amount of inclination. 

5. The oblateness of the earth also causes 
slight fluctuations of the radius (altitude) and 
period of tne satellite. 

All of the above effects are important in 
sustained satellite operations and tneir effects 
on orbits adds to the knowledge of the shape of 
the earth, however, such effects are of small 
magnitude in direct ascent rendezvous trajectories. 


As an example, in reference 6 the formula: 


2 
=> s 
A = Əw (1.637 x 107”) ———— cos i 


(R. + z.)” 

is given as the measure of angular motion of the 
nodes in radians per revolution of the satellite. 
For an altitude of 1000 n. miles and with a cir- 
cular orbit inclined at an arbitrary angle of 45°, 
the regression of the nodes is 4.38 x 16? radians 
for every revolution of the target. The ascent 
trajectories described in this thesis are of short 


duration compared with the period of the satellite 
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target and only a portion of this regression of 
the nodes takes place while the RV is in flight. 
For the 1000 n. mile target the ascent time varies 
from 3134 seconds to 843 seconds or 11% to 40% of 
the period of the target. Thus the regression of 
the nodes is from 13,000 to 47,000 feet at this 
geocentric radius. 

This example is given only to show the order 
of magnitude of the motion in this case. In the 
rendezvous situation the perturbations affect 
the RV as well as the target at the same time and 
in the same direction. The error is, then, of 
differential order when compared with the above. 

If these effects are expressed as functions 
of position they can be introduced into the inte- 
gration of the equations of motion for the 
rendezvous computations ,but they are not the sub- 
ject of the present investigation and can influence 
it only through higher order effects. 

To continue with the basic assumptions, 
booster burnout altitude 18 taken as sixty nautical 
miles for all cases. Without a detailed examination 
of the booster performance and maneuvering scheme 
1t cannot be determined whether this altitude is 


appropriate for a given trajectory,but some arbi- 
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trary altitude is needed for a common base from 
which to compare the coasting performance. This 
altitude is assumed in numerous references such 
as Eggleston in reference T. 

It is further assumed that the ephemeris of 
the target is known exactly. The technique of 
tracking an earth satellite and determining its 
ephemeris is fairly well formalized and nigh 
speed digital commune routines are available for 
` this purpose. The technique is described by Sandorff 
in reference 8. For this study, however, uncertainties 
of the satellite position are simply uncertainties in 
the origin of the final error coordinates. While this 
may be an important consideration, it is independent 
and can be introduced as necessary in a system study. 
Although the ability to predict launch time 
continually improves with experience in missile 
operations, there are still uncertainties in time 
of lift-off as well as error build-up during burn- 
ing. The time of flight or coasting time varies 
considerably with the particular trajectory 
selected (as seen below) giving some flexibility 
in matching time. A pre-launch computer keeps 


track of this time and, over some useful range 
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determined by other considerations, selects the 
exact trajectory to match time. In a similar 
mater, during the burning period performance is 
monitored and cut-off corrections computed to 
match the time. In very general terms, Duke, 
Goldberg and Pfeffer in reference 5 take notice 
of this type error but its implementation de- 
pends upon System Specifications. If an error 
in time alone is accepted and no compensation 
is made then the arrival error is easily com- 
puted from the satellite velocity. To get to 
the point of the study, it is assumed further 
that booster burnout is compelled to occur ex- 
actly on time. 

Uncertainties Considered 
The uncertainties considered in this thesis 
are those in the booster burnout velocity vector. 
These errors are described as: 
1. Position errors 
a. ahead (RANGE error) 
b. above (ALTITUDE error) 
c. to one side (TRACK error) 
2. Velocity errors | 


a. error in magnitude 
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b. error in direction 


(1) deflected above or below 
the ideal 


(2) deflected to one side of 
the ideal 


The exact meaning of "position errors" de- 
pends upon the type of guidance equipment utilized 
by the booster stage; the full implication of this 
point is carefully discussed in Chapter III. 
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CHAPTER III 


COMPUTATIONAL METHODS AND PROCEDURES 


Analytical Computation of Reference Trajectories 
The elements that describe the target orbit 

and the family of possible direct ascent ellipses 

tangent to it are defined with the nelp of Fig. 1 


below. 


V4 7 opodee velocity 


Ascent trajectory 







burnout velocity 






Surface 


Nominal limit Nef atmosphere 


Vs 
Circular satellite 
velocity 


Fig. 1. Definition of Elements Used in Computational 
Procedures. 
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From the figure it is seen that selecting 
a target altitude fixes R, for all ascent ellipses, 
since the rendezvous condition is that apogee 
occurs at target altitude. An additional variable 
is needed to specify a particular trajectory from 
the possible family of ellipses,and the range 
angle, d , best describes what is a long or short 
ascent. Although these two variables are selected 
for classification of results, the conventional 
two-body equations are more familiar in terms of 
the standard parameters e, a and 9. The desired 
values are computed by successive substitution in- 
to the numbered equations below. 

These equations are not involved, but they 
require careful numerical evaluation for accurate 
results. With dimensions in millions of feet, 
the IBM 650 carrying eight significant digits is 
adequate for computation directly in engineering 
units. The three physical constants for the 
numerical work are: 


G, (gravitational constant times mass 
of the earth) 


1.40768 x 1010 ft’/sec® Ref. 8 


R_ (radius of a sphere with volume 
equivalent to that of the earth) 


3437 .866 nautical miles Ref. 9 
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International nautical mile = 6076.1155 feet 








Ref. 10 
Altitude traversed in ascent 
Zu = 2. = 2 4 (3.1) 
Minimum eccentricity 
a ox a. | 
“min eR, TE, (5.2) 
Angle from perigee , 
O = (180° - db) (252) 
Eccentricity 
5 | 
cos9 + 1 l | cosQ + 1 гә /А, 
p — a lA A ед 
2] cos9+1+ — l | cog9+1+ =| |cos9+1+ — 
R R R 
1 1 1 
(5.4) 
Satellite velocity 
Gm 
Y-a S (3.5) 
S R; +2, 
Booster burnout velocity 
2 z 
V. = V LM. + l-e (3.6) 
1 8 R4 


Va se Y l-e 650) 
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Injection velocity increment 
V. = V. - Ç (3.8) 


Semi-major axis 


acu -а (2.9) 


Booster burnout attitude angle 


e sino 
€ cun Бе | (3104 


Time to apogee 


a - R a - R 
a 1 1 
T = a| am E - arccos (вә ) +e sin | arccos E Ih 
(@ 33.) 


This information, computed for 12 target alti- 
tudes from 200 to 20,000 nautical miles and range 
angles of 180, 90, 45 and 22.5 degrees permits 
selection of a smaller representative sample for 
detailed error analysis. Target altitudes of 400, 
1,000, 5,000 and 20,000 n. miles are used in the 


error analysis to follow. 


The Error Analysis 


A complete error analysis of even this basic 
set of trajectories by analytical methods, if 
possible at all, is certainly complex because 


several orbital elements are often changed by a 





21 
single error at booster burnout. For instance, 
the variation in semi-major axis can be found as 
a function of a change in Y , but before an error 
can be isolated, the corresponding change in ec- 
centricity and the shift of the line of apsides 
must be determined. Finally, arbitrary position 
as a function of time is trancendental. 

Braham and Skidmore in Hef. 11 give analytical 
expressions for changes in altitude, final velocity 
and final flight path angle as a function of certain 
initial condition errors. These expressions are 
ae. The one for final flight path angle 
requires evaluating six lengthy expressions to get 
the necessary partial derivatives, which in turn 
require evaluating four other equations equally 
as long to solve for some of the coefficients. 

Even if a complete set were available, evaluation 
in any quantity calls for a machine computer pro- 
gram. 

Such analytical solutions are interesting 
from the theoretical point of view, but it is 
felt that the most practical way to get a com- 
plete error analysis is by numerical integration 


of the equations of motion. 
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The motion of a body subject to an accelera- 


tion is described by the vector equations: 


V = fa dt + V (5.12) 
R= [Tae +F, (5.13) 
In coasting flight near the earth, the only 
acceleration is that due to the gravitational 
field, G, which is a function of R. Thus, for 
Step-by-step integration, the two simultaneous 


equations are actually solved in order as follows: 


iw Ei (3.14) 
R 

This acceleration is expressed in components in 

whatever coordinates V is given, and each com- 

ponent of V and R is computed as: 


x 


i AU T V 


= G kə (5-15) 


ü 


R, = V, At + Rio (3.16) 


i 1 


The resulting R is used to compute a new G; the 
resulting V is the V, for the next step, and the 
resulting R is also the Ro for the next step. This 
procedure is repeated up to apogee time. 

A complete program with this capability, de- 
veloped by Mr. Robert Brown of the M.I.T. Instru- 
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mentation Laboratory and adapted to this problem, 
15 used in the error analysis. With this program 
and the IBM 709 computer at the M.I.T. Computation 
Center, trajectories for four altitudes and four 
range angles are thoroughly investigated. 

The method is straightforward and accurate 
when used with & suitable integrating time step. 
The validity can be checked for each trajectory 
by integrating over the complete time of flight 
with no errors introduced and comparing the 
apogee conditions with the analytical results 
from the equations above. As an example of the 
accuracy attained the zero-error numerical in- 
tegration for ascent to a 5,000 n. mile altitude 
at a range angle of 90° using a 50 second time 
step agreed to within 80 feet ALTITUDE and 7 feet 
RANGE distance - this is at a geocentric radius 
Of over 51 million feet. 

The program uses spherical coordinates in 
the actual integration, but keeps track of the 
components in terms of an orthogonal frame called 
a platform frame. This frame ts introduced to 
simplify the presentation of the error components, 


and has a positive X axis in the RANGE direction 
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along the horizontal in the direction of flight, 
a positive Z axis upward along the vertical 
ALTITUDE direction, and a positive Y axis to the 
right looking in the direction of flight. 

A platform frame is specified because for 
inertially guided vehicles directions are re- 
ferred to a stabilized platform. The platform 
is non-rotating in earth-centered inertial co- 
ordinates (type 1) en its Z axis tracks the local 
vertical (type 2). The effect of position error 
on both types is shown i Pig. 2, 

The first type seems most likely. for this 
application, but it is used here for a more basic 
reason. With error effects available for the 
first type, they can be computed for the second 
type. Note from the figure that ALTITUDE errors 
alone are the' same in either system, but that 
witn RANGE or TRACK error, tne relative alignment 
of the G vector depends on the type System assumed. 
It is seen though, that a type 2 RANGE error con- 
sists of a type 1 RANGE error plus a velocity mis- 
alignment ($% ) which is easily computed from 
booster burnout radius and the RANGE error. 

Effects due to the actual change in height 
caused by a type 1 RANGE or TRACK error are negli- 


gible because, for reasonable errors, this change 
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in heignt is very small, and the corresponding 
position error derivatives are also very small. 
The same is not true for misalignment errors. 

The errors may be fractions of a milliradian, 

but this kind of error has a large end effect. 
Finally, for TRACK error, the motion due 

to the type 1 position error must be obtained 

by numerical integration, but for the type 2 

platform where the ALTITUDE component of velocity 

is aligned with the local vertical, the entire 

motion lies in the plane of the error trajectory. 

When a trajectory lies entirely in a plane and 

the inclination with respect to another plane 

is known, the lateral displacement is just a 

function of the geometry. This point is brought 

out in more detail in Chapter IV. 
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reference 
local 
vertical 







reference 
local 
vertical 


Comparison of two types of platform coordinate 
frames. Angle for reasonable position error 
greatly exaggerated. 





Xp 
Ret. local horizontal Ref. local horizontal 
Error in magnitude of In-plane misalignment 
booster burnout velocity. of booster burnout 
velocity. 
2 
reference ı 
ee | 
ҹеҹлиә. 
۷ 
$T d 
. V, cos Y 
4rack 
error 
Ye 
Track Error - Tp oe Lateral Misalignment 
system as used nerein. of booster burnout 


velocities 


Fig. 2. Platform Coordinates and Definitions of Booster 
Burnout Errors: 
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CHAPTER IV 


RESULIS 


summation of Velocities 
The summation of booster burnout velocity and 

injection velocity increment (ZV) for various alti- 
tudes and range angles is shown in Fig. 3. It must 
be borne in mind that the values represented here 
are tne theoretical minima; they do not reflect the 
additional velocity that must be provided to over- 
come drag and gravity losses during accelerations 
to Zu: 

| What is considered important about these plots 
is the effect of altitude and range angle on 2V. 
Note that while analysis always prescribes the 
Hohmann maneuver for minimum 27, this is not 
necessarily the primary consideration in a practical 
case. For instance, with a 400 n. mile orbit alti- 
tude the range angle is reduced to 55° before EV 
increases 1% over the Hohmann value; it is reduced 


to 18° before £V increases 5%. Similarly for a 
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1000 n. mile orbit altitude the range angle is 
down to 40° before €V increases by 5%. 
These examples are cited only to call at- 
tention to what may be one corner of a many 


sided compromise. 


Velocity Division f 
Fig. H is a plot of the separate components 


of Z V versus range angle for a number of orbit 
altitudes. Tne da is displayed on two 
separate plots for low and high altitudes to 
eliminate the clutter of several overlaps. 


V the booster burnout velocity, is the 


17 
velocity which must be imparted to the RV to 
send it coasting up to apogee. V5; the injection 
velocity increment, represents that propulsive 
capability that must be carried to the rendezvous 
point in order that satellite velocity be matched. 
From another point of view, Vo is the relative 
velocity between target and RV. If the ascent 
trajectory is perfect and time matched exactly, 
the injection thrusting can be initiated by elapsed 
time, but with error buildup during the ascent it 


seems most likely that tne terminal phase will be 





Velocity - ft/sec x 1077 


Velocity - ft/sec. x 10 7 
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— E Injection Velocity Increment 
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O 45 90 155 180 
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Fig. H. Effect of Altitude and Range Angle on Booster 


Burnout Velocity and Injection Velocity 
Increment. 
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closed-loop. Thus, this relative velocity may 
pë important during the acquisition of the esaet 
by tne RV. 
As before, for computation of propulsive re- 
quirements the drag and gravitation losses must 


be added to V.- Then, once an effective exhaust 


1 
velocity or specific impulse is known or assumed, 
quick preliminary calculations can be made con- 
cerning vehicle asdışa, These might consist of 
selecting an altitude and/or range angle vithin 
the capability of a certain system, or, in deter- 
mining vhat configuration and staging or vhat re- 
start capability is needed to accomplish a given 
mission. 
This information forms the second part of 


the propulsive picture, and another side of the 


overall compromise involved. 


Error Derivatives 

The results of the error analysis are re- 
duced to partial derivatives since it is shown 
that the variation of a final component with re- 
spect to a change in the initial value of any 
other component is linear over what is considered 


to be a generous range of variation. This direct 
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proportionality is verified for the 1000 mile 
orbit for speed errors up to 60 feet per 
second, for angular deviations of the velocity 
up to two milliradians and for position errors 
up to 12,000 feet in any direction. As shown 
for the errors in booster burnout velocity, 
plotted on Fig. 5, there is no tendency to de- 
part from linearity at 60 feet per second. Be- 
cause the functions being considered are smooth 
and continuous, they do not bend abruptly so it 
is considered that booster burnout velocity 
errors up to 100 feet per second and angular 
devlations involving components of that size 
could be entertained with confidence. It does 
not seem likely that such large errors are 
realistic, although what is tolerable certainly 
depends on what is said about the terminal 
system capabilities. In any case, the additional 
computer time required to investigate by this 
method and establish the absolute limit of appli- 
cability is not justifiable with the many otner 
items being studied with the computer. 

The types of booster burnout uncertainties 
described above are classifiable as in-plane and 


lateral errors. With three components each to 
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Fig. 5. Effect of Errors in tne Magnitude of the 


Booster Burnout Velocity on Range and 
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define the position and velocity vectors, six 
differentials are sufficient to describe any 
possible error in the booster burnout condition. 

To express the reference apogee or final value 
of each of tne six vector components in terms 
of each of the six initial components leads to 
36 partial derivatives. There is, however, 
clearly no way for the in-plane effects to in- 
duce any lateral perturbations and, as expected, 
the lateral errors introduce only infinitesimal 
effects on the final in-plane performance. With 
tnis independence, the number of error derivatives 
1s reduced to 20. 

While it seems more natural to describe tne 
initial velocity errors in terms of variations 
in magnitude and direction, the final velocity 
errors are given in terms of orthogonal components. 
This is done primarily for ease of combination of 
components as demonstrated later. Also, since botn 
the target and the reference apogee point have zero 
Y and Z components of velocity, the computed error 
components are the actual relative velocity in the 
ALTITUDE and TRACK directions as viewed from tne 
satellite. 


The X component of the error velocity, while 
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generally quite small compared to either 
the actual or relative velocities in the RANGE 
direction, is still. important. It is a part of 
the terminal problem, and the rendezvous is not 
complete until all components of relative 
velocity have been exactly matched. 

The derivatives relating the initial and 
final variations are made up of combinations 


of the following: 


1. In-plane.- Initial Position (x and z) 
and velocity (vi and Y) 


Final position (X and Z) and velocity 
(V, and V,). | 

2. Lateral - Initial Position (y) and velocity (í). 
Final position (Y) and velocity (Vy). 


The general form of the drivative as used for 
evaluating ascent errors is Зар. This 

does not preclude using the inverse. If, for 
instance, it is desired to define allowable final 
error first, then the permissible booster burnout 
deviations of a type could be found. However, 
there are, in general, many combinations of booster 
burnout errors that can combine to cause similar 


final errors. 
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The derivatives are plotted in Figs. 6 through 
15 with Pour angle as the primary variable 
and the target altitude as the step parameter. 
Each is explicitly defined on the respective 
plot with units in feet, feet per second and 
milliradians. 

1t is important to note that as range 
angle approaches zero, the eccentricity of the 
trajectory ellipse approaches unity daa the 
regular forms of the equations become very 
sensitive. Judging from a few of the cases 
for which the vertical ascent may be evaluated, 
it appears that OE SE functions in fact are 
quite unstable at low range angles. For this 
reason, none of the results should be extra- 
polated below the 22.5” range angle where the 
error computations were cut off. There seem 
to be practical reasons for avoiding very steep 
ascents such as the high relative velocities, 
highest possible XV, and need for rapid change 
of attitude in the terminal phase.: For such 
reasons, it is considered that die lack of 


coverage at these short angles is not important. 
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Independence of Derivatives 

A final important result is the independence of 
the in-plane derivatives. To investigate all possible 
combinations of booster burnout errors in varying 
magnitude involves endless computation, but the total 
error is shown in Table I to be made up of the sum of 
the errors due to the separate effects. Table I is 
compiled by introducing a random combination of errors 
into the 1000 n. mile, 90۶" range angle trajectory and 
checking the computed output for component addition 


from the separate derivatives. 
TABLE I 


Final Error Components 


Booster Burnout X Z Ys ə 
Error feet feet| ft/sec ft/sec 






(1) 420 ft/sec -44,2953 40.25 
(2) -2 mr -103,095 -10.03 
(5) 1 mile ahead -T,T11 -.75 





Sum of (1) + (2) 99,097 |-153.09 30.22 
98,929 |-155.10 29.89 
Sum of (1) + (3) -52,004 | 64,741 39.50 


Computed Results -51,994 | 64,753 39.48 


Computed Results 





Sum of 
(1) + (2) + (3) -155,099 | 105,145 |-164.01 29.47 
Computed Results | -155,080 104,964 |-164.04 29.09 
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Let tne assumed booster burnout errors be de- 


signated by $x, $z, SV. and é¥ for in-plane 


1 
errors and by ¿áy and éT for lateral errors. Nov 
1f the final errors are designated by AX, AZ, av, 
and AV, for in-plane errors and by AY and DN for 
lateral errors then a complete final error for a 
booster burnout error in, say, Speed, range and 


lateral alignment is: 


Ax = 9X. By A BZ $x 








ov, 1 Ox 
AZ = GA $v, + 2 $x 
AV = zə ÖV, 4 = Sx 
۵۷ = A oV. + s Sx 
AY = Sn 


o". c 
Avy "Se OA 


where the partial derivatives are read from tne 


graph for the appropriate altitude and range angle. 
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Table I demonstrates that this linearization 
introduces no more than 1% inaccuracy in the final 
computed values for a practical case. Special 
notation seems unnecessary since the procedure 
is so simple. Any combination is readily written 
if it jis remembered that each initial in-plane 
error causes four components of final error, two 
in position and two in velocity. The lateral 
errors, which are independent of the in-plane 
error, each cause two final errors, one in TRACK 


position and one in TRACK velocity. 


Lateral Position Error 

The lateral position error due to a TRACK 
error at booster burnout rates special mention 
because it illustrates an effect due entirely 
to the definition of position error used, that 
is, the type of platform assumed. The easiest 
place to view this effect is at a range angle 
of 90 degrees. It seems at first that for a 
lateral position error at booster burnout, the 
line of nodes between the plane of the actual 
ascent ellipse and the reference plane would 
lie 90 degrees away, or at apogee for this 
case. The integrated results shown in Fig. 14 


show this is not the case. 
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Fig. 16 below illustrates how this comes 
about with the central angle due to & reasonable 
position error greatly exaggerated. 


Reference plane 


Nore : Burnout osition and 
Velocity com ponents are 


In platform coordinates Y sin¥ (Y/R) 


ALTITUDE. component o$ V, 


M 
TRACK posit on error 1 


Xp into paper 


FA 


Fig. 16. Effect of Platform Definition on 
| Lateral Position Error, 
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Note that although no lateral velocity 
appears in the platform coordinates, with a 
platform that stays aligned with the reference 
vertical, a sidewise component of velocity, 


V. sim Y (WR, ); does show up with respect to 


1 
the reference orbital plane. 
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CHAPTER V 


CONCLUSIONS 


If tnis thesis has met its objective then 
the plots snowing tne various error derivatives 
Should clearly هوه‎ btw the primary con- 
clusions to be drawn concerning the effect of 
target altitude and range angle on these de- 
rivatives. 


To review the results in general: 


l. For tne low altitude rendezvous 
the summation of booster burnout 
velocity and injection velocity in- 
crement (ZV) is only slightly great- 
er than the Hohmann value until the 
range angle is well below 90°. The 
components of ZV follow this same 


pattern. 


2. The in-plane error derivatives 
for a 1000 n. mile target altitude 
are accurate for difforemtials of 
+ 60 ft/sec, + 2 milliradians and 
12,000 feet in any direction. 





D2 
2. The lateral error derivatives are 
independent of the in-plane error de- 


rivatives. 


4. Тһе in-plane partial derivatives 
are independent of each other and there- 
fore the total error components can be 


formed by the rule for total differentials. 


5. The propulsive requirements for 
the 20,000 n. mile orbit altitude are 
shown to be so great that multi-staging 
will certainly be necessary. 
At this altitude the sensitivity 
to error is so great that special guidance 


techniques will also be required. 
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